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We present a measurement of the top quark pair {tt) production cross section in pp collisions 
at \fs = 1.96 TeV using events with two charged leptons in the final state. This analysis utilizes 
an integrated luminosity of 224-243 pb~^ collected with the D0 detector at the Fermilab Tevatron 
Collider. We observe 13 events in the e+e~, e^ and At^/i~ channels with an expected background 
of 3.2 ± 0.7 events. For a top quark mass of 175 GeV, we measure a tt production cross section of 
Gtt = 8.6I2 ^(stat) ± l.l(syst) ± 0.6(lumi) pb, consistent with the standard model prediction. 

PACS numbers: 13.85.Lg, 13.85.Qk, 14.65. Ha 

The top quark was discovered in 1995 at the Fermi- Its observation completed the third quark weak isospin 
lab Tevatron Collider in pp collisions at ^/s = 1.8 TeV. doublet suggested by the absence of flavor changing neu- 
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tral current interactions [3| and measurement of the b 
quark weak isospin |^. By virtue of its large mass 
{nit = 178.0 ± 4.3 GeV Q), the top quark could de- 
cay into exotic particles, e.g. a charged Higgs boson [^. 
Such decays would lead to a measured tt production cross 
section {au) apparently dependent on the tt final state. 
It is therefore necessary to precisely measure att in all 
decay channels and compare it with the standard model 
prediction. The increased luminosity and higher collision 
energy of ^/s — 1.96 TeV at the Run II of Tevatron per- 
mit substantially more accurate measurement of au in 
all final states. 

In the SU (2) x U (1) electroweak model with one Higgs 
doublet 0] , each top quark of a ti pair is expected to de- 
cay approximately 99.8% of the time to a W boson and 
a b quark 0. Dilepton final states arise when both W 
bosons decay leptonically. These occur along with two 
energetic jets resulting from the hadronization of the b 
quarks and missing transverse energy (^t) from the high 
transverse momentum (px) neutrinos. In this Letter, we 
present a measurement of att with 224-243 pb^^ of pp 
collider data at y/s = 1.96 TeV collected with the up- 
graded D0 detector |^. We consider the e+e~, e^i and 
fj-^ fj-^ final states. The electrons and muons may orig- 
inate either directly from a W boson or indirectly from 
a.W ^ TV decay. The corresponding tt branching frac- 
tions [B) are 1.58%, 3.16%, and 1.57% for the e+e", 
e/x, and IJ.'^ i-J.^ channels, respectively. 

The D0 detector has a silicon microstrip tracker and 
a central fiber tracker located within a 2 T supercon- 
ducting solenoidal magnet |^. The surrounding liquid- 
argon/uranium calorimeter has a central cryostat cov- 
ering pseudo-rapidities up to 1.1 I^L and two end 
cryostats extending coverage to \ri\ « 4 [iflj. A muon sys- 
tem resides beyond the calorimetry, and consists of a 
layer of tracking detectors and scintillation trigger coun- 
ters before 1.8 T toroids, followed by two similar layers 
after the toroids. Luminosity is measured using plastic 
scintillator arrays located in front of the end cryostats. 
The trigger and data acquisition systems are designed to 
accommodate the high luminosities of Run II. The data 
used in this analysis were collected by requiring two lep- 
tons (e or n) in the hardware trigger and one or two 
leptons in the software triggers 0. 

To extract the tt signal, we select events with two high- 
Pt isolated leptons, large ^t, and at least two jets. We 
further improve the signal to background ratio by se- 
lecting events with kinematics compatible with events. 
To derive the cross section we determine the overall ef- 
ficiency e (including trigger, geometrical, and event se- 
lection efficiencies) for tt and the number of expected 
background events. We distinguish two categories of 
backgrounds: "physics" and "instrumental". Physics 
backgrounds are processes in which the charged leptons 
arise from electroweak boson decays and the origi- 
nates from high px neutrinos. This signature arises in 



Z/j* — > r+r" where the r leptons decay leptonically, 
and WW/WZ (diboson) production. Instrumental back- 
grounds are defined as events in which (a) a jet or a lepton 
within a jet fakes the isolated lepton signature, or (b) the 
originates from misreconstructed jet or lepton ener- 
gies or from noise in the calorimeter. 

The electrons used in the analysis are defined as clus- 
ters of calorimeter cells for which (a) the fraction of 
energy deposited in the electromagnetic section of the 
calorimeter has to be at least 90% of the total cluster en- 
ergy, (b) the energy is concentrated in a narrow cone and 
isolated from further calorimeter energy, (c) the shape of 
the shower is compatible with that of an electron, (d) 
the electron matches a charged track in the tracking sys- 
tem. In order to further remove backgrounds we use (e) a 
discriminant that selects prompt isolated electrons based 
on the tracking system and calorimeter information |l2l| . 
Electrons which fulfill criteria (a) to (e) are referred to as 
"tight" electrons. For background calculations we intro- 
duce "loose" electrons for which only (a) and (b) are re- 
quired. The muons considered in the analysis are defined 
as tracks reconstructed in the three layers of the muon 
system, with a matching track in the tracking system. 
The energy deposited in the calorimeter inside a hol- 
low cone around the muon must be less than 12% of the 
muon pt- To further remove background, the sum of the 
charged track momenta in a cone around the muon track 
has to be smaller than 12% of the muon pT- Muons that 
fulfill all these criteria are referred to as "tight" muons. 
For background calculations, we introduce "loose" muons 
for which the isolation criteria are relaxed. 

Jets are reconstructed with a fixed cone of radius 
ATZ = 0.5 01 and must be confirmed by the indepen- 
dent calorimeter trigger readout. Jet energy calibration 
is applied to the jets [lj| . The is equal in magnitude 
and opposite in direction to the vector sum of all signif- 
icant calorimeter cell transverse energies. It is corrected 
for the transverse momenta of all isolated muons, as well 
as for the corrections to the electron and jet energies. 

Event selections for each channel are optimized to min- 
imize the expected statistical uncertainty on the cross 
section. We select events with at least two jets with 
pip > 20 GeV and \y\ < 2.5 and two leptons with > 
15 GeV. Muons are accepted in the region \ri\ < 2.0, while 
electrons must be within |?7| < 1.1 or 1.5 < |7/| < 2.5. 
The two leptons are required to be of opposite signs in 
the e+e" and channels. 

A cut on is crucial to reduce the otherwise large 
Z/7* background. This background is particularly severe 
in the e"'"e~ and n'^ fJ-~ channels. Due to different reso- 
lutions in electron energies and muon momenta, the op- 
timization leads to different selections in the three chan- 
nels. In the efj, channel, we require > 25 GeV and 
A(j){^T,fJ') > 0.25, where A(j){^T,f^) is the azimuthal an- 
gle between the and the muon. The latter gives addi- 
tional rejection against Z/j* — > tt background in events 
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with two jets. In the e+e channel, we veto events with 
dielectron invariant mass 80 < M^e < 100 GeV and re- 
quire > 35 GeV i^T > 40 GeV) for M^e > 100 GeV 
(Mge < 80 GeV). In the fi~^iJ.~ channel, we accept events 
with > 35 GeV. This cut is tightened at low and high 
values of A(j>{^T , IJ-i) where /ii denotes the leading px 
muon. Events with A0(^t,Mi) > ^75° are removed. 

The final selection in the channel requires — 
pf^ + S(p^) > 140 GeV, where pfji denotes the pr 
of the leading lepton. This cut effectively rejects the 
largest backgrounds for this final state which arise from 
Z/"f* — > T^T^ and diboson production. The e^e^ anal- 
ysis uses a cut on sphericity S — 3(ei -I- e2)/2 > 0.15, 
where ei and £2 are the two leading eigenvalues of the 
normalized momentum tensor [l5| . This requirement re- 
jects events in which jets are produced in a planar geom- 
etry through gluon radiation. The final selection applied 
in the ijl^ IJ.~ channel further rejects the Z/^* —f 
background. We compute for each /x^/j," event the 
of a fit to the Z IJ''^l-i~ hypothesis given the mea- 
sured muon momenta and known resolutions. Selecting 
events with > 2 is more effective than selecting on the 
dimuon invariant mass for this channel. 

Signal acceptances and efhciencies are derived from a 
combination of Monte Carlo simulation (MC) and data. 
Top quark pair production is simulated using ALPGEN 
[iq with mt = 175 GeV. pythia is used for frag- 
mentation and decay. B hadron and t lepton decays 
are modeled via evtgen and tauola respec- 
tively. A full detector simulation using GEANT [20| is 
performed. Lepton trigger and identification efficiencies 
as well as lepton momentum resolutions are derived from 
Z/7* t^t' [Jl — e,/i) data. These per-lepton normal- 
ization factors and momentum smearings are applied to 
MC events to ensure the simulated samples provide an 
accurate description of the data. The jet reconstruction 
efficiency, jet energy resolution and -^t resolution in the 
MC are adjusted to their measured values in data. 

To calculate the expected number of events from 
physics backgrounds, we use Zj^* — s- T+r~ and dibo- 
son MC samples generated with pythia and ALPGEN, 
respectively. The Zj^* — > t'^t~ contribution is normal- 
ized to the cross section measured by D0 [ij]. For the 
diboson processes, diboson -I- 2 jets events are generated 
at leading order (LO) and are scaled by the ratio of the 
next-to-leading order to LO inclusive cross sections de- 
rived for diboson inclusive production |2^ . 

Instrumental backgrounds are determined from the 
data. Fake electrons can arise from jets comprised essen- 
tially of a leading 'k^ lr\ and an overlapping or conversion- 
produced track. We estimate this background by calcu- 
lating the fraction /e of loose electrons which appear as 
tight electrons in a control sample dominated by fake 
electrons. In the e+e^ channel the control sample con- 
sists of events that satisfied the trigger and have two 
loose electrons. In the e/i channel the events in the con- 



trol sample must satisfy the trigger and have one tight 
muon and one loose electron. Contributions from pro- 
cesses with real electrons {W — > ev and Z/7* e+e^) 
are suppressed by requiring < 10 GeV in both e^e~ 
and e/i channels and |Mee — Mz\ > 15 GeV in the 
e+e~ channel only. We also veto events in which both 
loose electrons have a matching track. We observe that 
/e measured in the e+e" and control samples agree 
within statistical errors. The predicted number of events 
with a fake electron in the final sample is obtained by 
multiplying the number of e~^e~ (e/i) events with one 
loose electron and one tight electron (muon) by /g. 

An isolated muon can be mimicked by a muon in a 
jet when the jet is not reconstructed. We measure the 
fraction of loose muons that satisfy the tight muon 
criteria in a control sample dominated by fake muons. 
In the channel the control sample is defined as 

events that have two loose muons. To suppress physics 
processes with real isolated muons the leading pT muon 
is required to fail the tight muon criteria. This cuts ef- 
ficiently Z/7* /i^/i^ events but also W ^ fif events 
where a second-leading muon might arise from a muon 
in a jet. The number of events with a fake muon con- 
tributing to the final sample is estimated by counting 
the number of events with one tight muon and a loose 
muon and multiplying it by In the e/x channel the 
contribution from events where both leptons are fake lep- 
tons is already accounted for by using fe . The remaining 
contribution from events with a real electron and a fake 
muon, is determined by combining fe and a fake rate 
obtained on a control sample that satisfies the e/i trigger. 

The processes Z/j* £~^£~ {£ = e,/t), while lacking 
high pt neutrinos, might have a significant amount of 
measured due to limited resolution. In the e"*"e~ 
channel, this background is estimated by measuring a 
misreconstruction rate on data and applying it to the 
simulation. We observe that the spectrum in e+e" 
events with 80 < Mee < 100 GeV agrees well with the 
spectrum observed in 7 + 2 jets candidate events. 
We obtain the misreconstruction rate in data as the 
ratio of the number of 7 -|- 2 jets events passing the 
selection divided by the number failing the selection. The 
misreconstruction rate is also consistent with Z/j* — > 
e+e^ + 2 jets simulation. This rate is multiplied by the 
number of events that fail the selections but pass 
all other selections. In the /i^/f~ channel, the expected 
contribution of Z/j* /i~''/i~ background in the final 
sample is derived from events simulated with ALPGEN. 
Good agreement is observed between the data and the 
simulation in the variables and A(j){^T , l^i) ■ This 
allows us to obtain the probability for a Z/j* — s- ^jl^ ^J.^ 
event to pass the selection from the simulation. The 
sample is normalized to the number of observed Z /^* — > 
/i+/i^ events in the data with 70 < M^^ < 110 GeV 
before the selection. 

The number of observed events and estimated physics 
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FIG. 1: Predicted and observed (a) number of events with 0, 1 and 2 or more jets with all other selections applied, (b) 
and (c) leading lepton pr in dilepton events after all selections. The Z/y* contribution includes e'^e~ , t^t~ e/i, and ii^ 
final states. The ti prediction is shown for a^t — 7 pb. 



and instrumental backgrounds in the dilepton + 2 jets 
sample, the integrated luminosities and the exB for the 
tt signal are given in Tablc^for each channel. We observe 
5, 8 and events in the e+e^, e/i and channels, 
respectively. We estimate the probability to observe > 5, 
> 8, and exactly events in the e~^e~, e/i, and fJ-^fJ-^ 
channels as 22%, 43%, and 5%, respectively, using the 
measured att and taking into account systematic uncer- 
tainties. By generating pseudo-experiments we estimate 
that 20% of the possible outcomes have lower likelihoods 
than that of our observation. The significance of the 
observed tt signal over the background is 3.8 standard 
deviations. 

To compute the cross section, we calculate in each 
channel the probability to observe the number of events 
seen in the data as a function of cr^j given the number 
of background events and the signal efficiencies. The 
combined cross section is the value of att that maxi- 
mizes the product of the likelihoods in the three channels. 
The resulting top quark pair production cross section at 
^/s = 1.96 TcV in dilepton final states is 

ati = 8.6t|?(stat) ± l.l(syst) ± 0.6(lumi) pb 

for mt = 175 GeV, within errors of the standard model 
theoretical prediction of 6.77 ± 0.42 pb and in agree- 
ment with the recent result in Ref. |2J| . We find a^t also 
consistent with measurements carried out in different fi- 
nal states ^2 ll^. The total systematic uncertainty is 
obtained by varying the background prediction and sig- 
nal efficiencies within their uncertainties and taking into 
account correlations. The dominant systematic uncer- 
tainties are given in Table |nj In addition a 6.5% sys- 
tematic uncertainty is assigned to the luminosity mea- 
surement li^. The top quark mass affects the signal ef- 
ficiency, resulting in a dependence of ati on rrit given by 
dati/dmt = —0.08 pb/GeV for mt in the range 160 GeV 
to 190 GeV. 



Figure n^a) shows that the observed number of events 
with 0, 1, and 2 or more jets, with all other selections ap- 
plied, is consistent with the prediction (assuming Uti = 7 
pb). Figure n^b) shows that the observed and predicted 
spectra after all selections agree well. Other kine- 
matic distributions in dilepton events are also well de- 
scribed by the sum of ti signal and background contribu- 
tions at various steps of the event selection. 

The leading lepton px spectrum in the ti dilepton fi- 
nal states has recently been studied by the CDF Col- 
laboration ^27^ and a mild excess has been observed at 
low transverse momenta. This is not confirmed by our 
data, as shown in Fig. ^c). To test agreement between 
data and the prediction, we generate pseudo-experiments 
from the predicted leading lepton pT spectrum and use 
our measured crti to normalize the ti signal. We find that 
31% of the pseudo-experiments are less consistent with 
the parent distribution than the data. We conclude that 
data agree well with the prediction. 

In summary, we have measured the top quark pair pro- 
duction cross section at -^/s = 1.96 TeV in e^e^, e/i and 
/i+/i~ final states to be Ott = 8.6t2;7(stat) ± l.l(syst) ± 
0.6(lumi) pb for rrit — 175 GeV, in agreement with the 
standard model prediction and with measurements in 
other final states. 
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TABLE L Expected signal (assuming mt = 175 GeV and 
(Jti = 7 pb) and background event yields for e^e~ , e/i, and 
fj,^ channels. Instrumental backgrounds include and 
fake lepton backgrounds. Total uncertainties are given. 



Channel 






Integrated luminosity (pb"'^) 


243 


228 224 


Physics backgrounds 
Instrumental backgrounds 
Total background 


0.3 ±0.1 
0.7 ±0.1 
0.9 ±0.1 


0.7 ±0.2 0.2 ±0.1 
0.2 ±0.1 l.llo3 
0.9 ±0.2 1.4 ±0.4 


exB (10~^) 
Expected signal 


1 q+0-2 


S.2tot 1.0 ±0.1 
5.lto l 1.6 ±0.2 


Total prediction 


2.8 ±0.3 


6.lt°oi 2.9 ±0.6 


Observed 


5 


8 


TABLE II: Summary of systematic uncertainties on au. 


Source 




Aatt (pb) 


Jet energy calibration 
Jet identification 
Muon identification 
Electron identification 
Trigger 
Other 




±0.8 -0.7 
±0.3 -0.6 
±0.5 -0.4 

± 0.3 
±0.3 -0.2 
±0.2 -0.3 


Total 




± 1.1 



gram. 
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